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Abstract. We show that the extremely high luminosity 
of the short-period recurrent nova T Pyx in quiescence can 
be understood if this system is a wind-driven supersoft 
x-ray source (SSS). In this scenario, a strong, radiation- 
induced wind is excited from the secondary star and accel- 
erates the binary evolution. The accretion rate is therefore 
much higher than in an ordinary cataclysmic binary at the 
same orbital period, as is the luminosity of the white dwarf 
primary. In the steady state, the enhanced luminosity is 
just sufficient to maintain the wind from the secondary. 
The accretion rate and luminosity predicted by the wind- 
driven model for T Pyx are in good agreement with the 
observational evidence. X-ray observations with Chandra 
or XMM may be able to confirm T Pyx's status as a SSS. 

T Pyx's lifetime in the wind-driven state is on the or- 
der of a million years. Its ultimate fate is not certain, but 
the system may very well end up destroying itself, cither 
via the complete evaporation of the secondary star, or in 
a Type la supernova if the white dwarf reaches the Chan- 
drasekhar limit. Thus either the primary, the secondary, or 
both may currently be committing assisted stellar suicide. 
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— novae, cataclysmic variables — stars: individual: T Pyx 
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1. Introduction 

The recurrent nova T Pyx is the most prolific of all recur- 
rent novae, with a mean outburst recurrence time of only 
22 years. But even away from outburst, T Pyx's observa- 
tional characteristics are remarkable. Most importantly, 
its quiescent bolomctric luminosity is L^ Q \ > 10 36 ergs s _1 
(Patterson et al. 1998). This is much higher than expected 
for a short-period cataclysmic variable (CV) like T Pyx, 
whose orbital period of P or b = 1.8 hr places it well below 
the period gap. 
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The implied disagreement with theoretical expecta- 
tions is severe. For reasonable estimates of T Pyx's system 
parameters, the standard model of CV evolution predicts 
mass transfer to be driven by gravitational radiation at a 
rate of M gr ~ 4 x 10 -11 Mq yr _1 . The corresponding ac- 
cretion luminosity is only L ~ 9 x 10 32 erg s _1 . Thus the 
gap between T Pyx's observed luminosity and the the- 
oretically predicted value spans at least three orders of 
magnitude. 

In order to account for T Pyx's extreme luminosity 
without having to postulate a similarly extreme mass- 
transfer rate for the system, it has been suggested that, 
even in quiescence, the system may be powered predom- 
inantly by nuclear burning on the surface of the white 
dwarf (WD), rather than by accretion (e.g. Webbink et al. 
1987; Patterson et al. 1998). However, even this can only 
reduce the luminosity gap by about an order of magnitude 
(c.f. Section |2^). Thus regardless of the actual mode of ra- 
diation energy release, T Pyx's extreme luminosity implies 
an accretion rate M acc > (4 - 40) x 10 -9 Mq yr _1 . This 
lower limit is in line with estimates based on theoretical 
nova models: in order to reproduce T Pyx's short outburst 
recurrence time, an accretion rate M acc > 10 -8 Mq yr _1 
appears to be required (Livio & Truran 1992; Prialnik & 
Kovetz 1995; Contini & Prialnik 1997). 

Clearly, T Pyx poses a fundamental challenge to our 
understanding of CV evolution. Put concisely, the ques- 
tion is this: why is T Pyx's accretion rate so abnormally 
high relative to "ordinary" CVs? The answer we propose 
in this Letter is that T Pyx has left the standard CV evolu- 
tionary track completely and is instead currently evolving 
as a wind-driven supersoft X-ray source. 



2. T Pyx as a Wind-driven Supersoft X-ray 
Source 

2.1. Motivation 

T Pyx's observational characteristics in quiescence are 
clearly extreme among CVs. However, as recently pointed 
out by Patterson et al. (1998) they are actually a good 
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match to the properties of the supersoft x-ray sources 
(SSSs). Despite this empirical connection, the standard 
model for binary SSSs - thermal-timescale mass trans- 
fer from a companion initially more massive than the 
WD - cannot work for T Pyx: the shortest orbital period 
found in the detailed study of this phase by Dcutschmann 
(1998) is about 6 hr, much longer than T Pyx's 1.8 hr. In 
fact, King et al. (2000) find that the two shortest-period 
systems among known binary SSSs - RX J0537.7— 7034 
(P orb ~ 3.5 hr) and SMC 13 (P orb ^ 4.1 hr) - are already 
difficult to account for within the thermal-timescale mass- 
transfer framework, even allowing for non-conservative 
evolution. 

We therefore propose that T Pyx is instead a mem- 
ber of an (even more) exotic class of SSSs - the wind- 
driven supersoft sources. The existence of this class has 
recently been suggested by van Teeseling & King (1998, 
hereafter VK98; see also King & van Teeseling 1998, here- 
after KV98). These authors show that if the secondary 
star in a close binary system like T Pyx is strongly irradi- 
ated by soft x-rays, a powerful wind may be driven from 
its surface. Under certain conditions, the mass and angu- 
lar momentum loss in this wind can dominate the binary 
evolution, which in turn can drive mass transfer at a rate 
comparable to the wind mass-loss rate. If the wind-driven 
mass-transfer rate is high enough to power the irradiating 
source, a self-sustaining, stable, wind-driven state is cre- 
ated. All that is required in this scenario is an event to 
trigger the wind-driven evolution. 

2.2. Theoretical Background 

We take as our starting point Equation 7 in VK98, which 
gives the expected mass-loss rate from a low-mass, near- 
main-sequence secondary star irradiated by intense soft 
x-ray emission from the vicinity of the primary as 

M w2 ~ -3 x l0- 7 ^(m 2J L 3 7?7 s ) 1/2 M^r" 1 . (1) 
an 

In this equation, m 2 — M 2 /Mq and r 2 — R 2 / Rq are the 
mass and radius of the secondary, an = a/(10 11 cm) is the 
binary separation and L37 = L/(10 37 erg s _1 ) is the lu- 
minosity of the irradiating source. The factor rj s measures 
the efficiency of the irradiating spectrum in producing the 
wind, and <p measures the area of the mass-losing regions 
on the secondary's surface, relative to ttR 2 . 

The irradiating luminosity is most conveniently pa- 
rameterized in terms of the luminosity expected for steady 
shell burning, L s b, which is 



L sh = 2.9 x 10 



37 



Ma. 



erg s 



(2) 



,10- 7 M Q yr- 1 , 

(Iben 1982). Defining L sb: 37 = L sh /{10 37 erg s" 1 ), we 
therefore write the irradiating luminosity as 



L37 = T]aL a b,37 = 2.9n a 



10- 7 M yr^ 1 



(3) 



and use the numerical factor rj a < 1 as a measure of 
the efficiency with which the accreted material is actu- 
ally converted into luminosity. Thus rj a — 1 for steady 
shell burning, whereas n a — L acc /L s b SsO.3 for accretion 
onto a massive WD in the absence of any nuclear process- 
ing (iacc = GMiM acc /Ri). Combining Equations [l] and || 
yields 



Mace = ~gM w2 



1.2 x lO-V^ 2 



VaVs (f^) 



2/3 (4) 

mi Mq yr \ 



Here, g — — M acc / M w2 is the dimensionless accretion rate 
(measured in units of the wind mass- loss rate), mi = 
Mi/Mq, and we have used Paczynski's (1971) approxi- 
mation for Rhila to recast r 2 ja\\ in terms of q. Here, 
Rl 2 denotes the volume-averaged Roche-lobe radius of the 
secondary. If a stationary, stable, wind-driven state exists, 
the accretion rate in it must be given by Equation ||. 

In order to determine g as a function of the system pa- 
rameters, one needs to calculate the mass-transfer rate due 
to Roche-lobe overflow (RLOF) in a semi-detached close 
binary system whose evolution is driven by a stellar wind 
from the secondary. This calculation is described in detail 
by VK98 and KV98. Briefly, such a system will quickly 
settle in a stationary state with RL2/RL2 — R2/R2, pro- 
vided such a state exists and is stable. The stationarity 
condition is sufficient to find g as 

(6/3 2 + 2a)-(5 + 3C)(l + a) 



.9 



(l + a)(5 + 3C-6a) 



(5) 



where £ is the effective mass-radius index of the secondary 
(describing its reaction upon mass loss) , and [3 2 is the spe- 
cific angular momentum of the escaping wind material, 
measured relative to the specific orbital angular momen- 
tum of the secondary. 

Equation [| gives the dimensionless RLOF mass- 
transfer rate in the stationary wind-driven state. It is 
identical to Equation 30 of VK98, except that we have 
retained the explicit dependence on (3 2 (VK98 set f3 2 — 1 
throughout). As shown by KV98, the stationary solution 
is stable for system parameters appropriate to T Pyx and 
M w2 oc L 1 / 2 oc Male - In deriving Equation ||, it has been 
assumed that all of the material in the stellar wind from 
the secondary escapes. Also, the effects of mass loss from 
the primary have been ignored. In reality, T Pyx may un- 
dergo significant mass ejection during its nova eruptions, 
and the long term average M w \ may be comparable to 
M w2 . However, we have verified that even strong, episodic 
mass loss from the primary has virtually no impact on our 
results for T Pyx in quiescence, unless the specific angular 
momentum carried away by the nova ejecta is extremely 
high (much higher than that of the primary) . In that case, 
the angular momentum loss associated with nova erup- 
tions would further accelerate (and could possibly even 
dominate) the binary evolution. We will return to this 
possibility in Section ^|. 
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We finally consider the orbital period derivative of such 
a wind-driven system. If this is entirely due to stationary 
mass loss and/or transfer, we can combine Kepler's 3rd 
law with Paczyhski's (1971) approximation for R^ 2 and 
differentiate logarithmically to obtain 

Porb _ 3C - 1 Mi _ (3C - 1)(1 + g) M w2 

Porb 2 Mi 2 M 2 1 ' 

(c.f. Equation 34 in VK98). However, the period deriva- 
tive measured in real systems may differ from this, since 
the timescale associated with such measurements is much 
shorter than the timescale on which the RLOF rate can 
adjust itself (trlof ~ ~ 10 * ^ where H I R * ~ 

10~ 4 is the scale height in the atmosphere of the sec- 
ondary near the inner Lagrangian point). Thus observa- 
tionally determined period changes could, for example, be 
due to fluctuations of M w2 or /3 2 on timescales shorter 
than trlof- 

2. 3. Application to T Pyx 

We are now ready to apply the wind-driven evolution sce- 
nario to T Pyx. For definiteness in our numerical esti- 
mates, we will adopt main-sequence-based values for the 
mass and radius of T Pyx's secondary: M 2 = 0.12 Mq, 
R 2 = 0.17 Rq. These follow from Kepler's law, the or- 
bital period, Patterson's (1998) power law approximation 
to the the M-dwarf mass-radius relationship of Clemens 
et al. (1998) and Paczyhski's (1971) approximation for 
Rl2 — R2 of a Roche-lobe filling secondary star. In addi- 
tion, we will use M\ — 1.2 Mq as an estimate of the white 
dwarf mass (Contini & Prialnik 1997). We therefore take 
the mass ratio in T Pyx to be q = M2/M1 — 0.1. We note 
from the outset that even this rather low value for q could 
be an overestimate, since T Pyx's wind-driven evolution 
may already have reduced M 2 below its main-sequence 
value. Our assumption that the secondary is still close to 
the main sequence amounts to saying that wind-driving 
has only just begun. 

We begin by noting that, for a low-mass secondary 
undergoing adiabatic mass loss, we may take C — — 1/3 
(VK98). Next, we consider two extreme estimates for the 
angular momentum loss parameter (3 2 . To obtain a lower 
limit, we note that the stellar wind material will carry 
away at least the specific angular momentum of the sec- 
ondary, in which case 0% = 1 (this is the case considered 
by VK98). On the other hand, the stellar wind material 
may extract angular momentum from the binary system 
by frictional processes. An upper limit to the amount of 
specific angular momentum that is likely to be extracted 
this way is given by the specific angular momentum of 
particles escaping through the outer Lagrangian points, 
which is j 2 ~ 1.65o 2 f2 (Sawada et al. 1984; Flannery & 
Ulrich 1977; Nariai 1975). Here is the angular velocity 
of the binary system, and, in our notation, the correspond- 
ing angular momentum loss parameter is (3 2 — 2. These 
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estimates, together with q — 0.1, yield g ~ 0.5 (/3 2 = 1) 
and g ~ 2 ((3 2 = 2). Substituting these values back into 
Equation ||, along with mi = 1.2, we obtain 

f 7 x 10-V^rfe M yr- 1 fa = 1 
M acc (T Pyx) ~ t (7) 
[ 1 x 10- Mq yr" 1 fa = 2 

for T Pyx if it has settled into a stationary, wind-driven 
state. The corresponding luminosity is given by Equa- 
tion |^ as 

C 2 x 10 3 erg s- 1 [3 2 = 1 

L(T Pyx) ~ \ (8) 
[ 3 x 10 37 (j) 2 r] s r]l erg s" 1 #2 = 2. 

Thus wind-driving can indeed account for T Pyx's extreme 
accretion rate and luminosity, provided that the various 
efficiency factors in Equations |?] and || are not too far from 
unity. 

VK98 have argued that r/ s ~ 1 in SSSs, since soft x- 
rays are absorbed well above the photosphere and should 
therefore be quite efficient at driving the wind. The fac- 
tor <j) is just the area of the mass-losing regions on the 
secondary divided by ttR 2 , so that mass loss from the 
entire front hemisphere would correspond to cf> = 2. We 
may therefore also expect (f> ~ 1, even if the secondary 
is partially shielded by an optically thick accretion disk. 
However, the most interesting parameter in this context 
is r\ a . As noted in Section [2.2| , energy release by accretion 
yields rj a < 0.3. The upper limit in this inequality corre- 
sponds to a Carbon-core WD of maximum mass (1.4 M Q ) 
and minimum radius (O.OO2P0; Hamada-Salpeter 1961). 
For a 1.2 Mq WD on the same mass-radius relation, we 
have 77 ~ 0. 1. Thus gravitational energy release alone can 
only meet the system's luminosity requirements if the WD 
is even more massive than we have assumed and (3 2 > 1. 

The alternative is that nuclear processing continues in 
T Pyx even in quiescence. This would imply that nuclear 
burning in T Pyx occurs both quasi-steadily (in quies- 
cence) and explosively (during outbursts), with the former 
taking place at a rate slightly below the accretion rate. 
From an empirical point of view, this does not seem un- 
reasonable: as noted by Patterson et al. (1998), some SSSs 
in the LMC and M31 are recurrent (see Kahabka 1995), 
and the galactic novae GQ Mus and V1974 remained lu- 
minous soft x-ray sources for several years after their nova 
eruptions (Shanlcy et al. 1995; Krautter et al. 1996). Sym- 
biotic stars (SySs) provide another interesting point of 
comparison. Several classical SySs, such as Z And, exhibit 
erupting behaviour even though though their quiescent lu- 
minosities (L ~ 10 3 i Q ) suggest nuclear processing as the 
dominant power source (Miirset, Nussbaumer, Schmid & 
Vogel 1991). In addition, the luminosities of symbiotic no- 
vae (a distinct class from the classical SySs) remain high 
for decades after their outbursts (Miirset & Nussbaumer 
1994). 
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From a theoretical point of view, the situation is some- 
what more difficult. The problem is that the rate at which 
steady nuclear processing can proceed is not, in principle, 
a tunable parameter. Based on 1-dimensional models, ex- 
plosive and steady processing are generally expected to be 
mutually exclusive regimes, whose dividing line is a func- 
tion of accretion rate and white dwarf mass (Iben 1982). 
Thus steady burning is not expected to occur for accre- 
tion rates less than M crit ~ 1.32 x 1Q- 7 M$$. The accre- 
tion rates predicted by our wind-driven model are below 
this line (albeit by only a factor of 2.5 for (3 = 2 and 
A%d = 1.2M Q ).On the other hand, nuclear processing 
on real WDs is unlikely to be well-described by spheri- 
cally symmetric models. For example, thermonuclear run- 
aways (TNRs) triggered by localized temperature or pres- 
sure fluctuations may not always spread and evolve into 
global TNRs (Shara 1982; Shankar & Arnett 1994). Quasi- 
steady quiescent burning in T Pyx could therefore conceiv- 
ably be the collective outcome of many successive local- 
ized TNRs. The individual mini-eruptions would not nec- 
essarily be obvious observationally, provided they are suf- 
ficiently small and frequent. These localised TNRs could 
process non-degenerate material more slowly than it ac- 
cretes, leading to the build-up of a more and more massive 
non-degenerate envelope and, eventually, to a global TNR. 

The preceding is a highly speculative scenario, and we 
do not mean to endorse it too strongly. We have outlined 
it mainly to provide a specific illustration of the general 
idea that (quasi-)steady and explosive nuclear processing 
might take place in the same object. This general idea is 
not new. It was first suggested by Webbink et al. (1986) 
and again by Patterson et al. (1998), both times with- 
out any specific model in mind. The main achievement of 
the wind-driving mechanism is that, given a high radia- 
tive efficiency, T Pyx's high mass accretion rate and lu- 
minosity can be accounted for self-consistently. We do, of 
course, acknowledge that the requirement r/ a ~ 1 implies 
a fair amount of fine-tuning, regardless of whether one 
invokes an extremely massive WD or quiescent nuclear 
burning. But some theoretical fine-tuning seems reason- 
able for a system like T Pyx, whose short orbital period, 
high luminosity and ability to produce nova eruptions are 
a unique combination among CVs and SSSs. This is not 
to say that the wind-driven evolutionary channel must be 
narrow: most other wind-driven systems may be charac- 
terized by somewhat higher accretion rates than T Pyx 
and may therefore be steady SSSs. Observationally, such 
steady, wind-driven SSSs may nevertheless be rare, since 
their evolutionary timescales would be even shorter than 
T Pyx's. 

We finally turn to the orbital period derivative. On 
substituting the values for g and M W 2 = —gM acc derived 
above into Equation 8 and inverting, we find that the 



expected timescale for period increase due to stationary 
wind-driven mass loss and mass transfer in T Pyx is 



Pa 



(T Pyx) 




1 



(9) 



fa = 2 . 



Patterson et al. (1998) derived P olb /P olh ~ 3 x 10 5 yr 
from the periodic dip in T Pyx's optical light curve, but 
more recent timings suggest a slower rate of change. In any 
case, as noted in the discussion following Equation |^, the 
observed P rb may not be a valid estimate of the stationary 
(long-term average) value. 



3. Discussion and Conclusions 

We have shown that the abnormally high luminosity of 
the short-period recurrent nova T Pyx can be explained 
if the system belongs to the class of wind-driven SSSs. 
T Pyx's accretion rate is then much higher than that of an 
ordinary CV, because the system's evolution is dominated 
by a radiation-induced wind from the secondary star. The 
wind-driven evolution is self-sustaining, because the wind- 
induced accretion rate is sufficient to power the luminosity 
that excites the wind. 

Our model predicts that T Pyx is an intrinsically 
bright soft x-ray source. T Pyx was observed serendipi- 
tously as part of the ROSAT all-sky survey, but not de- 
tected (Verbunt et al. 1997). This is not yet in conflict 
with a SSS model for T Pyx, because of the limited sensi- 
tivity of these observations and the relatively high column 
density towards the system (N-r — 2.1 x 10 21 cm" 2 ; Dickey 
& Lockman 1990[]). However, a sensitive, pointed obser- 
vation with Chandra or XMM could yield a detection and 
provide a direct estimate of the WD luminosity and tem- 
perature. For example, assuming a distance of 3 kpc (e.g. 
Patterson et al. 1998), the flux produced by a blackbody 
with L WD = 10 37 ergs/s and T WD = 2.4 x 10 5 K is be- 
low the ROSAT detection limit, but would provide around 
1000 counts in a 30 ksec observation with XMM. 

If T Pyx is, in fact, a wind-driven SSS, its current 
evolutionary timescale is r OV oi = M2/M2 = -Porb/-Porb ~ 
10 6 yrs. Clearly, this is an extremely short-lived evolution- 
ary state. So how will T Pyx's evolution actually end? We 
can think of at least three possible outcomes. First, the 
mass ratio of the system may be driven so low that irradi- 
ation and wind-driving cease to be effective. Second, the 
wind-driven evolution may proceed until the secondary 
completely evaporates, leaving T Pyx as an isolated WD. 
Third, if the total mass of the system is somewhat higher 
than we have assumed, the WD may reach the Chan- 
drasekhar limit and explode in a Type la supernova. 

The likelihood of the last outcome depends on whether 
the amount of material that is ejected in each of T Pyx's 
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recurrent nova outbursts is significantly less than that 
which is accreted during each quiescent interval. Theoret- 
ically, this is plausible (e.g. Prialnik & Kovetz 1995; Livio 
& Truran 1992). The question might also be addressed ob- 
servationally: the total mass accreted during quiescence is 
AM acc ~ MaccT r0 c ~ 1 - 20 x 10~ 7 M Q . The mass ejected 
during an outburst can be inferred by comparing the pre- 
and post-outburst orbital periods, although this is sensi- 
tive to the details of the ejection process (Livio 1991; Livio 
et al. 1991). Such a comparison could also distinguish be- 
tween wind-driven and nova-driven evolution scenarios. As 



noted in Section |2.2| , if the specific angular momentum of 
the mass ejected during T Pyx's nova eruptions is very 
high (jnova<^a 2 ^ >> ji), the angular momentum loss 
associated with nova eruptions could, in principle, itself 
drive T Pyx's high mass transfer rate. In the context of 
this model, the difference between pre- and post-outburst 
orbital periods must be relatively large. By contrast, the 
wind-driving model makes no such prediction. 

Given the strong possibility that the system will de- 
stroy itself on a short timescale, T Pyx represents an evo- 
lutionary channel by which short-period CVs may be re- 
moved from the general CV population. The existence of 
such a channel is desirable from a theoretical point of view. 
As most recently pointed out by Patterson (1998), the ob- 
served CV population shows a significant dearth of short- 
period systems relative to the predictions of standard evo- 
lutionary scenarios. Whether T Pyx's "assisted stellar sui- 
cide" represents a sufficiently wide channel in this context 
is an important question: do most CVs eventually enter 
a wind-driven phase or is T Pyx really a unique system? 
The answer clearly depends on how wind-driving is ac- 
tually triggered, an issue we have not addressed in the 
present work. One possibility is that triggering occurs dur- 
ing a period of residual nuclear burning in the aftermath 
of a nova eruption, though this would have to be limited 
to systems satisfying other constraints (e.g. low q and/or 
high M WD ). 
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